Flight  Mechanics  of  an 
Elastic  Synunetric  Missile 

by  Charles  H.  Murphy 
and  William  H.  Mermagen,  Sr. 


ARL-TR-2255 


April  2000 


20010604  084 


Approved  for  public  release;  distribution  is  unlimited. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other 
authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  thereof. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 
to  the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5066 


ARL-TR-2255  _ April  2000 


Flight  Mechanics  of  an  Elastic 
Symmetric  Missile 

Charles  H.  Murphy 

Weapons  and  Materials  Research  Directorate,  ARL 


William  H.  Mermagen,  Sr. 

Computational  and  Information  Sciences  Directorate,  ARL 


Approved  for  public  release;  distribution  is  unlimited. 


Abstract 


The  free-flight  motion  of  an  elastic  missile  is  approximated  with  three  bodies  connected  by 
two  massless  elastic  cantilever  beams.  If  the  mass  distribution  of  the  three  bodies  is  1-2-1,  the 
frequency  of  the  symmetric  oscillation  of  the  outer  bodies  is  within  5%  of  the  classical  frequency 
of  the  oscillation  of  a  free-free  beam.  A  second  combined  pitching  antisymmetric  flexing  motion 
can  occur  with  a  frequency  that  is  almost  twice  that  of  the  symmetric  flexing  motion.  As  the 
beam  stiffiiess  is  reduced,  the  symmetric  flexing  motion  frequency  approaches  the  rigid  body 
aerodynamic  zero-spin  frequency,  and  the  flight  zero-spin  aerodynamic  frequency  is 
considerably  reduced.  Moderate  beam  damping  can  cause  dynamic  instability  for  spins  greater 
than  the  zero-spin  aerodynamic  frequency.  Resonance  mode  amplification  can  occur  when  the 
spin  is  equal  to  the  zero-spin  aerodynamic  frequency,  but  more  importantly  it  can  occur  when  the 
spin  is  equal  to  the  two  elastic  flexing  frequencies.  Spin-yaw  lock-in  occurs  at  the  lower  elastic 
frequency. 
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1.  Introduction 


The  linear  flight  mechanics  of  spinning  projectiles  was  first  developed  after  World  War  I  and 
later  extended  after  World  War  n  [1-4].  These  results  were  reformulated  by  Flatus  [5]  in 
missile-fixed  coordinates  for  reentry  vehicles.  Thus,  the  linear  flight  mechanics  of  symmetric 
projectiles,  including  finned  missiles,  is  well  established.  Nonlinear  flight  mechanics,  as  well  as 
the  effect  of  hquid  payloads,  have  also  been  extensively  addressed  [4,  6,  7]. 

Recently,  with  the  development  of  long-finned  anti-armor  projectiles  and  long  flare- 
stabilized  reentry  vehicles,  the  effects  of  aero-elasticity  have  become  a  matter  of  possible 
concern  for  designers.  In  1992,  Flatus  [8]  developed  a  simple  theory  for  elastically  deforming 
missiles.  This  theory  showed  significant  decreases  in  pitching  frequency  when  the  smallest 
bending  frequency  was  less  than  four  times  the  pitching  frequency.  He  also  showed  that  linear 
beam  damping  could  actually  make  the  missile  dynamically  imstable  when  the  spin  was  greater 
than  the  pitching  frequency. 

Flatus  assumed  that  the  elastic  missile  oscillated  in  the  modal  waveform  of  a  nonspinning, 
free-free  beam.  For  uniform  mass  distribution,  the  waveform  is  symmetric  with  equal 
deflections  at  each  end.  The  aerodynamic  loads  on  a  finned  missile,  however,  are  quite  different 
at  the  nose  and  the  tail,  and  this  may  not  be  a  good  assumption.  Moreover,  it  would  be  desirable 
to  predict  the  resonant  response  induced  by  mass  or  aerodynamic  asymmetries  at  aerodynamic  or 
elastic  frequencies.  The  response  of  a  spin-stabilized  body  of  revolution  to  the  motion  of  an 
interior  elastic  beam-mounted  mass  was  studied,  and  beam  damping  destabilized  the  shell  [9, 
10].  Since  all  shells  have  spin  rates  that  exceed  their  pitching  frequencies,  this  is  in  agreement 
with  Flatus’  result. 

In  this  report,  the  beam  theory  [9,  10]  was  used  to  construct  a  simple  model  of  an  elastic 
projectile  acted  on  by  two  significant  aerodynamic  loads.  The  complicated  elastic  projectile  was 
replaced  by  three  rigid  projectile  components  connected  by  two  massless  cantilevered  beams. 
Next,  this  three-body  model  was  used  to  yield  good  predictions  of  the  frequency  of  the  motion  of 
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a  uniforai  free-free  beam  flying  in  a  vacuum.  With  the  inclusion  of  aerodynamic  terms,  the 
frequencies  and  damping  of  the  angular  motion  of  a  symmetric  elastic  missile  can  be  estimated 
from  a  5  X  5  determinant. 

Beam  damping  should  be  based  on  derivatives  in  a  missile-fixed  coordinate  system.  Thus, 
beam  damping  has  no  effect  on  the  trim  motion,  which  is  constant  in  these  coordinates. 
Resonance  frequencies,  however,  will  decrease  as  the  ratio  of  elastic  frequency  to  rigid-body 
pitching  frequency  decreases.  Beam  damping  can  adversely  affect  the  damping  of  transients. 
Finally,  the  effect  of  permanent  deformations  were  considered  (i.e.,  the  case  of  a  “bent”  missile). 
These  deformations  affect  the  size  of  resonant  trim  and  introduce  the  possibility  of  spin-pitch 
lock-in  [11]. 

2.  Linearized  Equations  of  Motion 


The  actual  projectile  is  approximated  by  a  homogenous  cylinder  of  fineness  ratio,  L,  and 
mass,  m,  divided  into  three  segments  with  fineness  ratios,  Ij,  and  masses,  mj.  A  fixed-plane 
coordinate  system,  which  pitches  and  yaws  with  the  center  segment  (component  1)  and  has  an 
origin  located  at  the  center  of  mass  (cm)  of  the  center  segment,  is  shown  in  Figures  1  and  2. 
(Fixed  plane  coordinates  roll  so  that  the  Z  axis  is  always  in  the  Xe-Ze  plane.)  Its  orientation  with 
respect  to  earth  fixed  axes  and  the  velocity  vector  {j3,a)  are  also  shown  in  the  figures. 

For  j  =  2  and  3,  let  {Xj,Yj,Z  j)  be  the  location  of  the  cm  of  component)  with  respect  to  the 
cm  of  component  1 ,  and  let  )  be  the  inclination  of  the  forward  body  (j  =  2)  and  of  the  rear 

body  (j  =  3),  as  shown  in  Figures  1  and  2. 

Next,  it  is  assumed  that  the  fore  and  aft  components  are  elastically  connected  to  the  center 
component  by  two  massless  cantilever  beams  of  lengths  a  and  b  and  stiffness  (EI)j.  The  moment 
and  force  exerted  on  these  components  by  these  beams  can  be  computed  from  simple  beam 
theory  [10].  For  motion  in  the  x-y  plane. 
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(1) 


and 

M,=b,jYj+b,jT^j,  (2) 

where  hu  =  -(12/a')(EI)2,  b,3  =  -(12/b')(EI)3,  b22  -  (6/a")(EI)2,  b23  =  - (6/b')(EI)3,  b32 
=  -(4/a)(EI)2,  and  b33  =  -(4/b)(EI)3. 


Similar  relations  apply  for  F^j  and  M^j . 

For  the  dynamics  analysis,  a  coordinate  system  with  its  origin  at  the  center  of  mass  of  the 
complete  projectile  is  much  more  convenient.  Throughout  this  report,  complex  quantities  are 
used  to  describe  the  lateral  motion  of  the  three  components.  Their  lateral  locations  are  specified 
by  the  dimensionless  complex  variables  Ej^  and  their  axial  location  by  the  dimensionless  variable 
Xj .  Both  the  lateral  and  axial  locations  are  with  respect  to  the  projectile  cm. 

Ej  - E,  =  {Yj  +  iZj )/ d;  E^  =  -(wjEj  +  "*3-^3 )/«■ ;  (3) 

Xj  ={Xjld)+x^;  and  x,  =-(/«2^2  +»t^X^)/md .  (4) 

The  complex  beam  forces  and  moments  exerted  on  the  forward  and  afl  components  are 


and 

9 

where  T,.  =  T  .  +  zT . . 

j  yj 


(5) 

(6) 
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The  velocity  vector  of  the  j -component  is  and  the  complex  transverse  velocity  is 

^yj  +  -  iXjQd,  (7) 

where  ^  =  /?  +  za ,  the  complex  angle  of  attack, 

Q  =  0  +  iify ,  the  complex  transverse  angular  velocity,  and 
V  =  v^j  is  the  magnitude  of  the  projectile’s  velocity. 

The  transverse  components  of  the  equations  of  motion  for  each  component’s  center  of  mass  can 
provide  differential  equations  for  the  Ej ’s, 

mj[Ejd*v{i-iQ)^.V4-ix,Qd\^Fj  +F.,  (8) 

where  Fj  is  the  transverse  aerodynamic  force  acting  on  the  j  component.  The  axial  beam 
forces,  F^j ,  are  determined  by  the  requirement  that  each  component  has  the  same  axial  velocity, 

=  (9) 

where  F^j  is  the  axial  aerodynamic  force  acting  on  the  j  component. 

The  angular  momentum  of  each  component,  Hj,  can  be  simply  expressed  in  fixed-plane 
coordinates  that  pitch  and  yaw  with  that  component, 

+  ^yj  )),  (10) 
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where  F,  =  0. 

After  differentiating  each  angular  momentum  vector,  three  differential  equations  can  be  written 
for  the  two  Tj  ’s,  and  Q  from  the  transverse  components, 

I,j(Q  +  ,fy)-ip/,j(Q  +  irj}=Mj  +Mj,  (11) 

where  My  is  the  transverse  aerodynamic  moment  acting  on  the  j  component. 

The  forces  exerted  on  the  center  component  by  the  other  components  are  the  negatives  of  the 
beam  forces,  and  the  moment  exerted  on  the  center  component  can  be  computed  fi-om  the 
negatives  of  the  beam  forces  and  moments. 

F,=-A-F,  (12) 

and 

M,=-M2-^3-i  [fe  -  +  fc  -  ^i)^3  +  ^2  te  -  Ej)+  F^3  {E3-  .E;)]  d .  (13) 

The  sum  of  equations  9  yields  the  usual  drag  equation, 

mV  =  'ZF„.  (14) 

y=i 

Equations  8  can  be  added  to  eliminate  the  beam  forces  to  provide  a  simple  relation  between 
0and(^, 


F{^-iQ)='LV,-4F,\ 

y=i 


(15) 
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Equations  8  is  now  multiplied  by  /(x^.  -x^)d ,  equations  9  multiplied  by  i{Ej  -E^)d ,  and  all 

are  added  to  the  sum  of  equation  1 1  to  eliminate  the  beam  forces  and  moments.  The  resulting 
differential  equation  for  Q  is  simplified  to 


ItQ-ipIxQ^  S  Vj  +  hj ry  -  ipixj  ] 

j=2 

+  I  [Mj+i[xjFj+F^jEjy[, 

J=i 


(16) 

where 


It  =It\  +It2  +It3  +  and 


I X  I x\  1x2  1x2  • 


3.  Linear  Aerodynamics 


The  linear  aerodynamic  forces  and  moments  acting  on  the  three  bodies  have  the  following 
forms  in  terms  of  the  angle  of  attack  at  each  body, 


=  -SxCoj, 


(17) 


I'yj  +  =  -^1  [  (p^ IV )%.  ^y  +  ^ jd I V 


(18) 


and 


M  j+iM^j=-ig^d 


"3y 


+  i{pd +c^jijd/V 


(19) 
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where  g,  =  pSV^  j 2 . 


The  Magnus  force  and  moment  coefficients  (c^jCjy)  are  usually  neglected  for  slowly 

spinning  firmed  missiles,  but  must  be  retained  for  spin-stabilized  bodies  of  revolution.  The  local 
angles  of  attack  at  the  three  components  are  determined  by  their  location,  inclinations,  and 
motion,  * 

ii~i-4xjdlv-r,+Ejdlv.  (20) 

The  transverse  aerodynamic  force  on  the  rigid  projectile  is  obtained  by  adding  the  three 
equations  18  and  using  equation  20  for  Tj  =  Ej  =  0. 

FyR  +  =  -gl  [(•=■  +  i(pd/vKk  +  and  (21) 

where  =  c] 

Similarly,  the  transverse  aerodynamic  moment  for  the  rigid  projectile  is  obtained  by  adding 
the  sum  of  times  each  of  equations  18  to  the  sum  of  equation  19  and  again  using 

equation  20. 

Myg^  +iM^R  =-igid  (c3  -\-i{pd/V)cl\  +  c^{^dlv)  ,  and  =c^  (22) 

In  reference  8,  the  relation  similar  to  equation  20  erroneously  uses  the  missile-fixed  derivative  of  Ej,  and  not  the 
derivative  in  a  nonspinning  coordinate  systenL 
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where  -  <^3  -  X!  +  ^jC\j  ]  and  =  C4  =  ^  [^4^  +  yc j  (c2j  -  c^j )-  (xy )  ^Cj ;]. 

1  1 

From  equations  14  and  17,  the  familiar  form  of  the  drag  equation  is  written 

mV  =  -g,C^,  (23) 

where  —  c^j  +  0^,2  +  . 

For  the  aerodynamic  force,  equations  17-19  can  be  inserted  in  equation  14,  and  the  small 
effects  of  the  Magnus  force,  the  damping  force,  and  the  center  body  lift  on  the  trajectory  can  be 
neglected  {c[j  =  =  c„  =  0) . 

mv[i-iQ)=-g^  -4c//F)-c^(r  -4c//F)j,  (24) 

where 

Equations  8  and  11  for  j  =  2.3  and  equations  15  and  16  are  six  differential  equations  in  six 
variables — Q,  E2,T2,E^,  and  If  the  aerodynamic  forces  and  moments  of 

equations  17-19  are  inserted  and  Q  is  eliminated  by  equation  24,  a  set  of  five  equations  in  five 
variables  results,  as  in 

E  [KJ.  +  (S„„  +  ipSl  )z„  +  (r„„  +  ipT*^  )z„  ]  =  0 ;  and  n  =1 , 2,  3,  4,  5,  (25) 

where  {Zi,Z2,Z2.z^,z^)  =  {^,E2,T2,E2,T3)  ■ 
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and  are  primarily  dynamics  terms  involving  various  masses  and  moments  of  inertia. 
The  aerodynamic  contributions  to  these  terms  are  quite  small  and  will  be  neglected. 
contains  combinations  of  beam  elastic  terms  and  aerodynamic  terms,  while  and  are 
combinations  of  beam  damping  terms  and  aerodynamic  terms.  In  equation  25, 


S..=L*ig,‘l/V)S^.md  (26) 

TL=fL-<-(g,d/yfL- 


Twenty  of  the  dynamics  terms,  twelve  of  the  elastic  terms,  and  eight  of  the  beam  damping  terms 
are  nonzero;  these  forty  terms  are  given  in  Table  1.  The  forty-six  nonzero  aerodynamic  terms 
are  given  in  Table  2. 

A  general  solution  of  equation  25  can  be  expressed  as  the  sum  of  ten  modal  functions  of  the 
form 


Zn,=  ZmicC^*  and  A  =  Ale  ,k=  1  to  10.  (27) 

When  a  function  of  this  type  is  substituted  into  equation  25,  the  final  result  can  be  written  in  the 
form  of  five  linear  homogenous  equations  in  the  constants  .  These  equations  are  specified  by 

a  5  X  5  matrix  (Unm)  which  is  a  function  of  A. 

tv-=0  andn.1,2,3,4,5.  (28) 

m=l 


where  u„„  =  {s„„  +  ipSl„ )+  +  ipT^„ . 
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The  ten  values  of  Ak  are  the  roots  of 


det  Unm  0.  (29) 

Forz,^  56  0 ,  the  corresponding  values  of  the  Zmk’s  for  each  of  the  ten  modes  can  be  determined 
from  the  solution  of  a  fourth-order  inhomogeneous  linear  system, 

4 

Z  ^nmymk  =  aud  n  =  1 , 2,  3,  4,  (30) 

m-\ 


where  “  U(n+ 1  )(m+ 1 )? 

ymk  ^  Z(ni4-i)k/ and 

Wn  =-U(n+l)l. 

Equation  29  can  be  solved  in  general  by  a  simple  trial-and-error  Newton’s  quadratic  formula. 
The  first  of  the  five  modal  equations  is 


I,A^+{Sjj-ipI,)A  +  Tjj 


ipTjj 


=[. 


7  7 

m2(x2  -Xj )d  A 


(lf2A  -  ipl^2  U  -  -  T,3  -  ipT^3  T2 


+ 


Xi)d  A  Sj^A  Tj^ 


ys 


^is-^  '^15 


\y4- 


(31) 


(The  k  subscript  has  been  omitted  from  equation  31  for  simplicity).  The  rigid  projectile 
frequencies  and  damping  rates  can  be  computed  from  equation  31  when  the  yj’s  are  zero.  Very 
good  approximations  [4]  for  the  frequencies  are 


^kR  = 


Ph 


<kr*«,Tu 


+  (l±pIj4I,as)pIJ2I„s^ii  =1,2,  (32) 
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where  is  the  rigid  projectile  zero-spin  frequency.  Exact  relations  for  the  damping 

rates  are 

K  =-(^,,4  -PJ,]-  (33) 

The  magnitude  of  ^2*  specifies  the  relative  size  of  the  angular  motion  of  the  fore  component 
with  respect  to  the  projectile  angular  motion  for  mode  k,  while  the  magnitude  of  specifies 
the  angular  motion  of  the  rear  component  for  that  mode.  For  example,  the  argument  of 

determines  whether  the  fins  of  an  elastic  missile  are  more  or  less  effective  than  those  on  a  rigid 
projectile.  Usually,  the  phase  angle  is  near  zero  and  they  are  less  effective. 

4.  Three-Component  Motion  in  a  Vacuum 

A  nonspinning  projectile  with  identical  components  2  and  3  and  identical  connecting  beams 
is  now  considered.  For  this  case,  all  parameters  for  components  2  and  3  are  equal,  except  for 
Xj  =  -Xj  and  bjj  =  -^23  •  For  motion  in  a  vacuum,  gj  =  0  and  the  aerodynamic  terms  vanish. 
The  five  equations  of  motion  become 


- »l2X2{z2k  - Z4k)^^  - +  25* =  0.  (34) 

-m2X2A^Zii,+[m2A^-bi2(l  +  r^)]z2k  -(^22/^)^3/t -^nVAk  (35) 

/,2^'zu  +  )z2i  +  (-  -^,2^'  +  ^32  K  +  ^  >  (36) 

m2X2A^Zik-bi2r^Z2k+lm2A^ -b^2(}  +  r„)]p4,,+{b22/d)z5k=0,  (37) 
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and 


-biidr,z^^  -b^dfy  +  r,)z,^+(-I„A^  +ft„)z5j  =0, 


(38) 


where  r^=m^lm^. 

Equation  34  indicates  that  zero  is  a  double  root  (Ai  =  A2  =  0).  If  modal  functions  are  considered 
for  which  Z2k  =  Z4k  and  zsk  =  -zsk ,  zik  is  zero,  and  equations  37  and  38  are  identical  to  equations 
35  and  36.  Thus,  four  more  roots  are  the  solutions  of  a  quadratic  in  A^,  as  in 

[m^A^  (39) 

For  small  7,2 , 

4(3  4  =  ±  /  ^2)mEI /{m^m^a^  j,  and  g  =  ±  .  (40) 

The  third  and  fourth  roots  are  particularly  interesting.  Planar  motion  occurs  when  Z23  =  Z24 .  The 
sum  of  the  modal  functions  for  Ej  is 


E2  =E^  =  2z23  cos  0)  t. 


(41) 


where  S  =  I{A^  }  and  l{a  +  ib]  =  b . 

This  in-phase  symmetric  motion  of  the  fore  and  aft  components  is  quite  similar  to  that  for  the 
classical  jfree-fi'ee  beam  vibration.  Geradin  and  Rixen  [12]  give  the  theoretical  fi'equency  for  two 
lower  free-free  vibrations  as 


OJ^=Xl^Ell(mL^d^}, 


(42) 
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where  =4.730,^2  =7.853. 


If  the  length  of  the  beam  is  assured  to  be  the  distance  between  the  cm’s  of  the  center  component 
and  the  forward  component. 


and 


a  =  ;c,4  =  (i/2Xl +  '•.)/(! +  2'-.). 

(43) 

\ld>  =  (4.566X/r„  (l  +  )V(l  +  2r„  )'  . 

(44) 

6?,  /d)  is  plotted  vs.  the  mass  ratio  in  Figure  3,  and  the  three-body  approximation  is  quite 

good  for  fore  and  aft  components  which  are  one  half  the  mass  of  the  center  component  or  are 
each  25%  of  the  total  mass.  In  the  numerical  examples  of  this  report,  each  end  component  is 
specified  to  have  one  quarter  of  the  total  mass. 

The  remaining  four  roots  can  be  obtained  by  assuming  an  antisymmetric  oscillation  of  the 
end  bodies  (i.e.,  zak  =  -Z4k  and  zsk  =  zsk).  Equation  34  becomes 


^ik  = 


2m2X2d  Z2k-^Iti^zk 


I  It, 


(45) 


and  equations  37  and  38  reduce  once  again  to  equations  35  and  36.  The  last  four  roots  are 
solutions  of  the  following  equation  in  A  : 

+b^^\  =0,  (46) 


where  a,  =  1  -  Im^  {x^  Y  I  If  and  21, ^  / 1,  • 
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A5,,=±i73El/(m,a,a^),  and  =±2i 


(47) 


The  ninth  and  tenth  roots  are  approximately  equal  to  the  seventh  and  eight  roots.  These  roots 
correspond  to  very  high  frequencies  and  will  not  be  considered  in  the  remainder  of  this  report. 

If  /,  is  assumed  to  be  the  value  for  a  slender  uniform  cylinder  of  length  L,  /,  is  ml}  Jxi. 
The  ratio  of  to  A^  ,^  can  be  derived  from  equations  40, 43,  and  47,  as  in 


^5,6  =  4.4  (l  +  2r„  )/Vl  +  2ri  .  (48) 

Roots  five  and  six  represent  a  pitching  motion  of  the  center  component,  combined  with  an 
opposing  oscillation  of  the  fore  and  aft  components.  For  =  0.5 ,  this  pitching  antisymmetric 

motion  is  80%  faster  than  the  free-free  in-phase  oscillation  that  was  considered  earlier.  The 
second  free-free  oscillation  is  antisymmetric,  but  is  180%  greater  than  the  first.  Thus,  the  three- 
body  theory  qualitatively  describes  this  oscillation,  but  predicts  its  frequency  very  poorly. 

5.  Flare-Stabilized  Rod 

Flatus  [8]  first  applied  his  theory  to  a  5-ft-long  rod  stabilized  by  a  very  light  1.5-ft-long  flare 
with  a  base  diameter  of  0.8  ft.  Since  he  neglected  the  very  small  normal  force  associated  with 
the  nose,  this  missile  has  a  very  simple  normal  force  distribution  of 

dC 

— —  =  -1.138{x  +  1.25)  and  -3.12  < x<-L25 .  (49) 

ax 

The  displacement  oscillations  of  a  free-free  uniform  beam  can  be  described  by  the  following 
sum  of  eigenftmctions: 


Mx,t)  =  (x) exp{i<j) and  —L{2<x< Lj 2 . 


(50) 


The  <5^  ’s  are  constants  determined  by  initial  conditions.  The  first  two  eigenfrequencies,  o,^ , 
were  given  by  equation  42,  and  the  corresponding  eigenfunctions  are 

^Xjc)  =  cosh(>v;c/Z)-T/^cos(\j[:/L),and^2(A:)  =  sinh(X2;r/Z,)-T/^sin()^x/Z,), 

[pi) 

where  r]f.  =  7.527, 35.975 . 

The  other  pairs  of  eigenfunctions  for  K  greater  than  2  have  the  same  form  as  equation  51,  with 
larger  values  of 

The  linear  motion  of  an  elastic  missile  with  distributed  aerodynamic  loads  can  be  expressed 
as  a  sum  of  periodic  functions.  Their  frequencies  are  related  to  the  aerodynamic  frequencies  of 
the  rigid  projectile  and  the  free  oscillation  frequencies  of  the  elastic  projectile.  Flatus  assumed 
that  each  of  the  periodic  components  had  the  form 

Although  the  first  three  eigenfunctions  were  used  to  compute  the  motion  of  the  flare- 
stabilized  rod,  the  second  and  third  eigenfunctions  had  very  little  effect  on  the  low  frequency 
motion.  The  first  eigenfunction  gave  the  following  expression  for  the  lowest  four  frequencies: 

'iM-piJ,  -Iii][fef +4M] = (53) 

The  ratio  of  the  displacement  at  position  x  and  the  angle  of  attack  for  the  k-th  frequency  is 


Z (^)  exp(4t)  and  l{4 }  =  .  (52) 
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-1 


(54) 


sMU = +/.//.] 

Table  3  provides  the  appropriate  parameters,  including  Flatus’s  integrals  .  The 

parameters  for  the  three-body  theory  are  given  for  w,  /2 . 

The  ratio  of  the  three-body  elastic  frequency,  d>,  to  the  rigid-body  zero-spin  pitch 
frequency,  ,  is  a  good  measure  of  the  effect  of  elasticity  on  the  projectile’s  aerodynamic 
frequencies  and  damping  rates.  This  ratio,  G)j(o^  ,  is  denoted  by  cr  and  is  usually  greater  than  5. 
When  it  is  less  than  5,  a  large  effect  on  the  aerodynamic  motion  is  expected.  By  varying  El, 
i^R  computed  as  a  function  of  cr .  The  result  is  plotted  as  Figure  4a  and  compared 

with  the  three-body  theory.  In  Figure  4b,  the  two  deflection  magnitudes,  of  the  three- 

body  theory  are  also  plotted  against  a,  and  they  are  compared  with  the  corresponding  deflection 
magmtude,  [<5j(3Z,/8)]/^ ,  of  the  Flatus  [8]  theory.  (The  three-body  deflections  are  given  by  y,, 
and  in  equation  30.) 

These  figures  show  similar  behavior  for  the  two  theories.  The  three-body  theory  predicts  a 
much  slower  decrease  in  ,  with  decreasing  <7  and  a  displacement  of  the  aft  component  that  is 

twice  that  of  the  forward  component.  Forex  =  2 ,  the  flare  angle  of  attack  is  reduced  by  40%, 
thereby  reducing  the  frequency  by  23%. 

6.  Finned  Missile 

The  implications  of  this  aeroelastic  theory  can  best  be  seen  by  considering  a  particular 
hypersonic  finned  projectile  with  a  fineness  ratio  of  20  flying  at  6,000  ft/s.  This  projectile  has  a 
body  diameter  of  0.35  ft,  and  its  forward  and  rear  segments  are  each  25%  of  the  total  length;  this 
is  slightly  less  than  the  optimum  length  of  27%.  The  various  aerodynamic  coefficients  are 
computed  from  an  assumed  force  distribution  (see  Appendix  A).  All  necessary  parameters  are 
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listed  in  Table  4.  For  zero  spin,  the  rigid  body  frequencycy^  is  8.52  Hz,  and  the  rigid  body 
damping  rate  is  -2.4  s'*.  For  cr  =  8,  the  three-body  symmetric  beam  frequency  co  is  69.0  Hz, 
and  the  elastic  projectile  frequency  is  8.06  Hz  (i.e.,  6.5%  less  than  the  rigid  body  frequency). 

The  equations  derived  in  section  3  can  now  be  used  to  calculate  as  a  function  of 

C7  between  1.0  and  10.0.  The  result  is  plotted  in  Figure  5,  while  the  magnitudes  of  and 
Y-^l^dxe  plotted  in  Figure  6.  For  cr  =  4,  the  effective  angle  of  attack  of  the  fins  has  been 

reduced  by  30%,  and  the  adverse  angle  of  the  nose  increased  by  35%,  thereby  reducing  the 
projectile  frequency  by  45%.  The  large  nose  deflection  with  its  increased  destabilizing  moment 
has  a  large  effect  on  the  projectile  frequency,  which  did  not  appear  for  the  flare-stabilized 
missile.  In  Figure  7,  the  elastic  frequencies  and  ^sI^r  plotted  against  cr.  The 

frequency  of  the  antisymmetric  pitching  oscillation,  ^5 ,  is  seen  to  be  approximately  twice  ^3 ,  the 
symmetric  oscillation. 

7.  Beam  Damping 

Beam  damping  is  a  very  complicated  process  that  will  be  approximated  by  a  simple  linear 
proportionality  with  the  displacement  velocity.  This  velocity  must  be  computed  in  a  coordinate 
system  spinning  with  the  beam.  Thus,  the  damping  force  exerted  by  each  beam  is  assumed  to 
have  the  form: 

(4  +'4L  -Ey‘*]=-2k,m,d>\E, -E,  -ip(E^  -£,)], and 
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where  =  ^3{EI%/[m^m^b^)  and^  =  p. 

The  definition  of  kj  was  selected  so  that  kj  =  1  corresponds  to  the  critical  damping  of  a  free- free 
beam.  For  example,  the  forward  beam  damping  contribution  to  the  first  mode  is 

-  [a,  4-  i(^^  -  p)\  [(w,  +  +  ^374,  ]/  m, . 

Thus,  beam  damping  is  a  strong  function  of  spin  and  is  essentially  zero  when  spin  is  equal  to 
the  modal  frequency.  The  beam  damping  terms  are  given  in  Table  1(c).  In  Figure  8,  the  finned 
missile’s  fast  mode  damping  rate.  A, ,  is  plotted  vs.  spin  for  cr  =  5  and  k2  =  ks  =  0,  0.05,  and 
0.10.  Beam  damping  increases  projectile  damping  when  spin  is  less  than  the  modal  frequency 
and  decreases  projectile  damping  when  it  is  greater  than  modal  frequency.  Figure  9  shows  A2, 
the  second  projectile  damping  rate,  for  the  same  parameters.  For  this  mode,  the  coning  motion  is 
opposite  to  the  spin,  the  frequency  is  approximately  the  negative  of  the  first  frequency,  and  beam 
damping  always  increases  the  damping  of  this  mode.  Finally,  Figure  10  shows  A,  vs.  <r  for  spin 

of  2.3(0j^,  and  k2  =  ks  =  0,  0.05,  and  0.10.  For  low  values  of  <7 ,  beam  damping  can  actually 
cause  projectile  instability.  Indeed,  Chadwick  and  Murphy  [9,  10]  show  that  the  beam  damping 
of  an  internally  mounted  beam  can  destabilize  a  spinning  shell  (Appendix  B). 

8.  Bent  Missile 


If  the  beam  deflection  is  large  enough,  the  maximum  stress  exceeds  the  yield  limit,  and  the 
beams  are  permanently  deformed  or  bent.  All  future  elastic  motion  will  be  about  this  new  bent 
shape.  Thus,  the  component  location  variables  Ej  and  are  no  longer  completely  elastic,  but 

consist  of  an  elastic  part  and  a  permanently  bent  part  that  rotates  with  the  projectile, 

Ej=Ej+Ej,e‘*  (56) 
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and 


r  =f  +r  e'^ 

1  -  i  -1- 1  jge  , 


(57) 


•  A  /K 

where  (j)  =  p  and  £'y,ry  vary  elastically. 

Both  fore  and  aft  components  will  have  trim  force  and  trim  moment  terms  associated  with 
and  Ejs  in  equation  20.  If  the  fins  are  deformed  relative  to  the  aft  component,  additional  fin 
force  and  moment  terms  must  be  added  for  the  rear  component,  as  in 

(58) 

and 

(a/„  +  M„  =  ig,dc„T„e* ,  (59) 


where  and  are  measmes  of  the  deformation. 

The  equations  of  motion  derived  in  section  3  should  be  modified  by  inserting  the  new 
variables  of  equations  56  and  57  in  the  dynamics  and  aerodynamics  terms,  adding  the  fin  trim 

terms  to  the  aerodynamic  force  and  moment,  and  replacing  Ej,Vj  with  Ej,tj  in  the  elastic 
terms.  The  new  set  of  five  inhomogeneous  differential  equations  contains  terms  multiplied  by 
exp(i^) ,  as  in 

Y [kJ.  +  +  ipSL K  + 1  +  ipTL K]= ^»exp(r^)  and  n  =  1, 2, 3, 4, 5,  (60) 

m=l 

where  t„  =  [p^  -ip)t,^  +  P^t*^  +git„A  +ip(SidlvY^ 

Kd  YodnAYnA  dcfincd  in  Table  5 .  Equation  60  can  be  solved  by  the  assumptions  of  constant 
spin  and  a  solution  with  the  following  form: 
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(zj ,  Zj  .^3  r  -24 » ^5 )  -  (^r  ’  ■^27- » ^27- ,  -^ST-  >  ^37- )  Qx^{ipt)  =  (5, ,  ,  ^3 ,  ^ 4 ,  )  exp(z>/) .  (61) 


The  Sj’s  are  solutions  of  the  following  linear  inhomogeneous  system  of  equations: 


’  n  =  1,  2,  3,  4,  and  5,  (62) 

m=\ 

where  =  -p^ (r^  +  Si, )+  ip{s„  +7;,)+  T„ . 

A  simple  bent  version  of  the  finned  missile  is  considered.  All  of  the  bent  parameters  are 
taken  to  be  zero,  except  for  £33  —  -.063  and  £3^  =  .020.  In  Figure  11,  the  zero-spin  trim 

Koi ^  divided  by  its  zero-spin  rigid  body  value, | ,  is  plotted  against  a.  For  cr  =  4 ,  the 

flexible  missile  has  a  trim  angle  twice  that  for  a  rigid  missile,  and  this  angle  increases  rapidly 
with  decreasing  cr .  The  magnitude  of  the  trim  motion  is  smaller  than  that  for  zero  spin  motion, 
except  in  the  vicinity  of  the  five  resonant  spins.  For  cr  =  5 ,  the  maximum  values  of  the  trim 
motion  parameters  for  spin  near  the  three  lower-resonant  frequencies  are  given  in  Table  6,  and 
Sj/sio  vs.  spin  is  plotted  in  Figures  12(a)-12(c)  for  j  =  1,  3,  5. 

For  spin  near  0.7 and  4.9®^,  the  forward  and  aft  components  move  in  phase  with  each 
other  while  they  are  out  of  phase  for  spin  near  9.7cOjf .  This  out-of-phase  motion  is  similar  to  the 
fifth  and  sixth  modes  for  vacuum  flight.  For  resonant  spins  of  O.To^j  and  9.7®^ ,  the  aft 
component  moves  to  reduce  the  effective  angle  of  attack  of  the  fins,  but  it  actually  increases  the 
effective  angle  of  attack  for  the  intermediate  spin  of  4.9® .  The  motion  of  the  two  end  bodies 
does  dominate  the  two  elastic  beam  resonant  fi'equencies. 
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9.  Quadratic  Roll  Equation 


Since  the  trim  motion  is  large  in  a  very  small  region  near  resonant  spin,  the  trim  motion  of 
rigid  spinning  finned  missiles  is  usually  a  minor  concern  for  designers.  Some  missiles,  however, 
can  “lock  in”  at  resonant  spin  for  a  range  of  initial  conditions.  A  rigid  missile  with  mass 
asymmetry  and  a  trim  aerodynamic  force  can  lock  in  at  resonant  spin  [11].  Quadratic  terms  in 
the  missile’s  roll  equation  cause  this;  therefore,  the  quadratic  version  of  roll  equations  must  be 
derived  for  the  three  components. 

Each  component’s  roll  equation  is  expressed  with  respect  to  that  component’s  fixed  plane 
X-axis, 

hjPj=M,j+g,c,jd.  (63) 

If  linearity  is  assumed,  then 

Pl=P2=P,=‘P,  (6^) 

(65) 

and 

^xP  =  ,  (66) 

where  c,  =c,,+c,,+c,,=C,^{p- p,,)dlV . 

is  the  linear  steady-state  spin  produced  by  either  intentional  cant  of  the  fins  or  unintentional 
damage  to  the  fins. 
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The  exact  relations  for  p2  and  ps  depend  on  the  requirement  that  they  should  be  equal  to  the 
corresponding  body’s  fixed-plane  X-component  of  the  central  body’s  angular  velocity. 

Pi  =  p't-r]  l2]+0  =  P+RiQ-pTil2f,) ,  (67) 


where  =jr^.|, 

+  ib}  =  a ,  and 

z  is  the  complex  conjugate  of  z. 


Equation  65  can  be  replaced  by  a  more  exact  expression  for  the  roll  moment  exerted  on  the 
central  body  by  the  two  beams;  this  replacement  is 


-I^,(E,-E,}l  +  F,{E,-E,)d]. 


Equations  63,  67,  and  68  can  be  combined  to  give  the  nonlinear  version  of  equation  66,  resulting 
in  equation  69. 


k -[hirl  +hiYl)lAp  =  SxCtd  +  R^2r2  +^2^3]-I^2{E2  -E,)d  +  FXE2  -EiVj-Q,  (69) 

where  Q  =i.hjRiQ- pFj /Tfj]  . 

2 


Equation  69  is  the  nonlinear  spin  equation  and  must  be  combined  with  equation  60  to  describe 
the  motion  of  a  bent  missile.  Any  study  of  the  steady-state  solutions  would  require  the  numerical 
integration  of  equations  60  and  69.  The  linear  spin  equation  had  a  constant  steady-state  solution, 
and  constant  steady-state  solutions  for  the  nonlinear  spin  equation  are  possible.  Indeed,  the 
nonlinear  spin  equation  for  a  rigid  bent  missile  has  only  constant  steady-state  solutions  [11].  For 
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constant  spin,  the  trim  motion  given  by  equations  61  and  62  describe  the  steady-state  transverse 
motion  of  an  elastic  missile.  If  the  equations  for  trim  motion  are  inserted  in  equation  69  and  p  is 
taken  to  be  zero,  the  appropriate  equation  for  equilibrium  spins  can  be  written. 


S^C„(p-p„)d^|V  =  -{MX.  (70) 

where  =  r\^ 2t^2t  ^ 2T^iT\~ V 

The  trim  values  of  Ej  andF^  are  defined  in  Table  5(b).  Two  functions  of  spin,  /j(p)and  /jCp), 
are  defined  as 


Irf ,  ^if2{j>hp-Ps.-  (71) 

Equation  70  states  that  equilibrium  spins  are  spins  for  which  /,  intersects  .  /,  {p)  is  a 
function  of  spin,  and  except  for  near  resonance  spin,  it  is  quite  small;  equilibrium  spin  must  be 
near  the  design  .  Near  resonance,  however,  /,  {p^  is  large.  Although  could  be  quite  far 

fi-om  resonance  spin,  a  near  resonance  spin  could  satisfy  equation  70  and  be  an  equilibrium  spin. 
If  this  equilibrium  spin  is  a  stable  equilibrium,  it  would  be  a  possible  constant  steady-state  spin 
and  therefore  produce  spin-yaw  lock-in.  Numerical  integration  of  the  differential  equations  of 
motion,  equations  60  and  69,  would  determine  the  stability  of  the  equilibrium  spins. 

In  Figure  13,  /i(p)  for  a  nonlinear  roll  moment  is  computed  fi'om  equation  71  for  the  bent 
missile,  and  the  result  is  plotted  against  plo)j^  .  The  line  /2(p)  for  the  linear  roll  moment  is  also 
plotted  for  p^  —IcOj^.  The  three  intersections  of  the  two  curves  are  equilibrium  spins.  The 
effect  of  other  values  of  p^^  can  be  determined  by  changing  the  p/cOr  intercept  of  the  /2(p) 
line.  This  process  shows  that  0.7®^  could  be  an  equilibrium  spin  when  p^^  is  between  0.76>;j 
and  3.1(»jj ,  and  similar  observations  apply  for  4.9®^ . 
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The  equations  of  motion  for  the  bent  missile  case  were  integrated  for  and 

pQ  -4.9(0^  ,  5.2  Q)j^ .  The  resulting  spin  history  is  plotted  against  time  in  Figure  14.  For 
greater  than  5.2<2)^ ,  the  spin  goes  to  a  constant  value  near  the  linear  prediction  of  p^^ .  For  p^ 
between  0  and  4.96)j^ ,  it  goes  to  a  value  of  4.7(o^  near  the  lowest  elastic  frequency.  Figure  14 
shows  that  the  boundary  value  of  p^ ,  which  separates  these  different  behaviors,  is  approximately 
5.05fu^  .  Thus,  spin  lock-in  can  occur  for  this  elastic  frequency. 

10.  Summary 

(1)  The  aero-elastic  motion  of  a  flexible  projectile  is  approximated  by  three  bodies 
connected  by  two  massless  elastic  beams. 

(2)  This  three-component  projectile  theory  gives  excellent  estimates  of  the  first  elastic 
frequency  of  a  free-free  uniform  beam. 

(3)  Two  types  of  beam  motion  are  simulated — ^the  classical  free-free  symmetric  motion  and 
an  out-of-phase  antisymmetric  flexing  motion  combined  with  pitching  motion. 

(4)  When  the  elastic  frequency  is  less  than  eight  times  the  aerodynamic  rigid  body  pitch 
frequency,  the  fins  move  to  reduce  the  stabilizing  moment  significantly,  and  the  pitch 
frequency  can  be  significantly  reduced. 

(5)  Moderate  beam  damping  can  cause  dynamic  instability  when  the  spin  is  much  greater 
than  the  pitch  frequency. 

(6)  The  trim  motion  of  a  bent  missile  shows  large  resonant  responses  when  the  spin  is  near 
the  lower  three  modal  frequencies. 

(7)  The  quadratic  roll  equation  shows  that  spin-yaw  lock-in  can  occur  for  the  lower  elastic 
frequency. 
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Figure  1.  X-Y  Coordinates  for  Three-Component  Projectile  (Every  Variable  Is 
Shown  as  Positive,  Except 


Figure!.  X-Z  Coordinates  for  Three-Component  Projectile  (Every  Variable  Is 
Shown  as  Positive,  Except  F^ ). 


25 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2 

Figure  14.  vs.  Time  for  Bent  Missile  With  cr  =  5  (i?,,  =  7©^;  pp 

=  4.9c0r;  5.2(«)r). 
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Table  1.  Nonerodynamic  Matrix  Elements 


(a)  Dynamics  Tenns 

It 

R4,  =  -m^Xjd 

R,2  =  -Xj)d^ 

R44  =  nijd 

1 

II 

^5/  “  ”“^55  “  ^t3 

Rj2 

Sn  =-h 

1 

11 

Sn  =1.2 

R2I  =-^2^2^^ 

Sis  =1.3 

R22  =  fn2d 

Si.  =  <3  =-1.2 

R3J  =  —R33  =  1^2 

Sl.=-Sl3=-I.3 

(b)  Beam  Elastic  Terms 

^22  =-b,2{mi+m2)dlm, 

T42  =-b,2m2d/m, 

i 

11 

T44  =-b22{m,-¥m2)dlm, 

^24  =-b,2m2dJmj 

T45  =  “^25 

T32  =b22{mj+m2)d/m, 

T52  =b23m2dlm, 

•>C) 

II 

T54  =b23{f»,  +m2)djm, 

T34=b22m2dlm, 

^55  “  ^33 

Beam  Damping  Terms 

S22  =  -T22  =  2k2m2G)2  {m,  +ni2)dlmj 

S42  =  —7^2  =  2k2m2G)2in2  djm, 

S24  =  -T24  =  2k 211226)2  m^dlnij 

S44  --T44  =  2k2m3G)3{m,  +m^djm2 
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Table  2.  Aerodynamic  Matrix  Elements 


^51  ~  ^55  ~  ^33^ 


u: 

1 

> 

1 

II 

^33  ■"  (p42  St2^12 

^12  ““(^32  ■*■^2^/2)^ 

'-i 

*N 

1 

11 

^13  ~  ^42  ^2^22  ~  St^I2 

S41  =  {C23  -  X3CJ3  -  m2X2  Q„  / 33l) 

=“(^55 

f 

11 

^/5  ~  (^43  “*■  ^3^23  ""  ^/^/3  )^ 

S^2=m2X2  0,2/m 

S2,  ={C22  -^2^12  -^2^2  Cia!^) 

S44  =  {m- 3312)0,2 /m 

S22  ={m -3312)0,2  Im 

S4S={^3X3/^)Ci3-<^23 

^23  “  (^2-^2  /^)^/2  ^22 

^51  ~  ~i^43  “  ^3^33  ~S t3^La 

^2.  = -m2  0,2!  331 

^53  ~  '^St3^I2^ 

S22  =3312X20,2/331 

S54  = 

^32  “  ”^32^ 

»n 

1 

II 
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Table  2.  Aerodynamic  Matrix  Elements  (continued) 


s„,„ 

where  g,  =  ijmd^ ;  andg,3  = 

-  I^^/md'^ . 

T* 

nm 

^35 

^3  “  (S^;c^l2  ^32  ^2^12  V 

^51  ~“(^jr34a  ~^3i)d 

^15  “  (^;c^l3  ^33  ^3^13 

^53  “  Sx3^l2^ 

=-{Sx2C,a  -4V 

'^SS  ~  (S^jc3^13  ”^^33)^ 

^3  —  (^a:2^12  ^32 

where  gy,=lj md^;  g^2  =  hi! >  and  g^-^  =  43 / md'^. 
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Table  3.  Flare-Stabilized  Rod  Parameters 


m  =  2.40  slug 

a  =  b  =  1.875ft 

m2  =  m3  =  0.60  slug 

p  =  0.002  slug/ft^ 

Ix=  0.0132  slug-ft^ 

V=  18,000  ft/s 

1x2  =  1x3  =  0.0033  slug-ft^ 

i 

II 

0 

II 

0 

It  =  5.00  slug-ft^ 

C3  =  -5 

It2  =  Ic  =  0.0781  slug-ft^ 

C33  =  -0.32 

L  =  5ft 

all  other  c ’s  are  zero 

X2  =  -X3  =  2.34  cal 

Ii  =  1.176  X  lOMb/ft 

(EI)2=  (EI)3 

I4  =  2.550  X  10^  Ib-ft 

d  =  0.80  ft 

l6  =  4.609  X  lO'’  Ib/ft 

drod  =  0.20  ft 

Ml  =  18.36  slug 
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Table  4.  Finned  Missile  Parameters 


m  =  3.50  slug 

Ct53  =  0.18 

nij  =  mj  =  0.875  slug 

Cl  =9.4 

=  0.054  slug  -  ft^ 

C3  =  -34.4 

1x2  =  1x3  =  0.0135slug-ft^ 

C4  =-750 

It  =14.20  slug-ft^ 

C12  =  2.3 

hi  -  Its  =  0-22  slug-ft^ 

c,3=7.1 

^22  ~  I  ^ 

(EI),  =  (e4 

C23  =  27 

t/  =  0.35ft 

C32=5.4 

a  =  6  =  2.63ft 

^33  “  ““3.6 

p  =  0.002  slug/ft^ 

C42  =  -30 

V -6,000  ft /s 

C43  =  *“10 

Cd  =  0.40 

C43 = -10 

C]32  =  0.15 

all  other  c-  are  zero 
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Tables.  Bent  Missile  Terms 


(a)  Values  of  C 


t\D  ~  ^\2^2B  '*'■^13^25  •^14-^35  ^35 


^2D  ~^22^2B 


11 

■^33  ^2B 

^4D 

^44^35 

ho  ” 

^55  ^35 

Ha  = 

-%2^1B 

“^3^25  - 

^14-^35  “^5^35 

~  [^33r3A/  +  ^3^13^3^]^ 

hA  ^ 

Jr  r 

'^23^  23 

^25(^35 

r„) 

t3A  = 

-C^2^2B^ 

11 

=  -T43^ 

^2B  ~  ^45 

(Fsb  +  r3p) 

hA  ““ 

’^^33^33 

+r,„)4 

4=0 


3D  ~  ^33^  2B 


423=0 


*  -  c*  p 
5D  *“  ^55^  3B 


^22^2B  ~  ^23^2B  ^24^33  ^25^33 


^^42^23  ^43^23  ^44^33  ^45^33 


(b)  Trim  values  of  Ej^Tj 


EzT  *^2  ■*■  ^23 


E^T  “■  *^4  ^33 


^27  *“  'S'3  +  T2S 


*^5  ^35 


Table  6.  sj/sio  for  Resonant  Spins,  cr=  5 


p 

0.74®  jj 

4.87®j, 

9.72c0j, 

mag  (si/sio) 

1.2 

HHDSHHH 

mag  (S2/S10) 

4.9 

22.1 

■KSlIHHii 

mag  (S3/S10) 

2.2 

8.8 

19.5 

mag  (S4/S10) 

4.7 

16.2 

93.3 

mag  (S5/S10) 

7.3 

20.6 

arg  (S]/sio) 

-87.3° 

106.0° 

0 

0 

0 

1 

arg  (S2/S10) 

93.3° 

102.6° 

-103.8° 

arg  (S3/S10) 

92.7° 

102.4° 

-103.6° 

arg  (S4/S10) 

91.8° 

102.4° 

75.6° 

arg  (S5/S10) 

-87.8° 

-77.5°  1 

-104.6° 
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Appendix  A: 

Aerodynamic  Properties  of 
Finned  Missile 
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Our  finned  missile  is  assumed  to  have  uniform  mass  distribution  over  a  20-cal.-long  cylinder 
with  the  density  of  aluminum.  The  external  configuration  has  a  fineness  ratio  of  22  due  a  very 
light  conical  windshield  of  1  cal.  and  a  1-cal.  extension  of  the  four  fins  past  the  cylinder  base 
(Figure  A-1).  The  centers  of  mass  of  the  three  components  are  located  at  xj,  and  these  locations 
can  be  computed  in  terms  of  the  lengths  of  the  cylindrical  segments. 


Figure  A-1.  Sketch  of  a  Hypersonic  Missile. 


x,=(i3-4)/2  ;x2=(i:-4)/2  -,x,^-{L-L,)l2.  (A-1) 

In  Wood  and  Murphy,'  the  slender  distribution  of  normal  force  coefficient  is  shown  to  be: 


cICn 

2 

dA 

/If 

djAd 

9d 

_  dx  _ 

S 

I 

_  ^  dx 

V 

V  . 

-10<x<ll, 


(A-2) 


) 


where  A  is  the  body  crossectional  area  and  Uj  is  angle  of  attack  at  Xj . 


*  Wood,  R,  M.,  and  C.  H.  Murphy.  “Aerodynamic  Derivatives  of  Both  Steady  and  Nonsteady  Motion  of  Slender 
Bodies.”  Journal  of  the  Aeronautical  Sciences,  vol.  22,  pp.  870-871,  December  1955;  see  also  BRL-MR-880, 
U,S.  Army  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  MD,  April  1955. 
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For  the  flight  dynamics  of  this  report,  d  =  d  and  the  9  term  have  no  effect. 


dC^ 

2 

dA 

/y _ 

dj  d  ^ 

dx 

body  ^ 

"  dx 

L 

V  1 

(A-3) 


This  expression  usually  underestimates  by  30%  due  to  its  prediction  of  no  normal  force 
acting  on  the  constant  area  cylinder.  Therefore,  we  add  a  cylinder  normal  force  term. 

=  0.6  \i.  -  (x  -  X.  \ddlV )  ]  ;  x  <  10 .  (A-4) 

cylinder 


dx 


The  contribution  of  the  fins  can  be  estimated  by  the  “strip”  theory  of  Boltz  and  Nicolaides,^ 
which  is  based  on  the  two-dimensional  linear  supersonic  theory  for  the  two  fins  in  the  x-z  plane. 


dC^ 

dx 


Jins 


=  2{aj^d^/s){y^  -0.5)[a^.  -(x 


Xj){6cd/V% 


(A-5) 


where 


Yf  = -(l3  +  3x)/4  .— 7<x<-5, 

=  2  -ll<x<-7 


Ul  -  2 


Vm^-1 


-1  -1 


M  =  Mach  number. 


[dCtialdx)  is  the  sum  of  three  coefficients  of  Uj  in  equations  A-3  through  A-5.  (C^y)^.  can  be 
obtained  by  integrating  the  sum  of  equations  A-4  through  A-6  over  the  i-th  component. 


-5L 


dx 


(A-6) 


^  Boltz,  R.  E.,  and  J.  D.  Nicolaides.  “A  Method  of  Determining  Some  Aerodynamic  Coefficients  From  Supersonic 
Free  Flight  Tests  ofa  Rolling  Missile.”  Journal  of  the  Aeronautical  Sciences,  vol.  17,  pp.  609-629,  October  1950. 
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(Cn),  =  J  ^  &=c„ft+c^(ftrf/r). 


(A-7) 


(C,),=  |r^]&  =  c„ft+c,,(|,rf/F).  (A-8) 

-iiL  ^  - 

(C„);  can  now  be  obtained  by  multiplying  the  sum  of  equations  A-3  through  A-5  by  (r- 

and  integrating  the  result  over  the  j-th  component.  In  Table  4,  these  relations  were  used  to 
compute  Cy  for  L2  =  L3  =  L/2. 

Finally,  the  roll  damping  moment,  was  computed  from  {dCi^ldx)j.^  using  the  local  angle 
of  attack  due  to  spin. 

c„  (pdiv)  =  -4{a,d' Is)  -  0.sjpi^^  +  0.5)/2r}fc .  (A-9) 

-11 
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Appendix  B: 

Internal  Beam-Supported  Mass 
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Chadwick'  and  Murphy^  studied  the  effect  of  beam  damping  on  the  dynamic  stability  of  a 
spinning  projectile  with  a  beam-mounted  symmetric  mass  and  aerodynamic  damping.  These 
references  showed  that  beam  damping  could  cause  dynamic  instability  when  the  beam’s  elastic 
frequency  is  near  the  projectile’s  fast  pitch  frequency.  The  effect  of  an  internal  cantilever  beam- 
supported  asymmetric  mass  on  the  motion  of  a  spinning  projectile,  however,  was  first  studied  by 
Reed^  for  no  aerodynamic  or  elastic  damping.  He  showed  that  dynamic  instability  would  occur 
when  the  elastic  frequency  was  near  the  projectile’s  slow  pitch  frequency  and  that  large  resonant 
motion  would  occur  for  spins  near  the  cantilever  beam  elastic  frequency. 

The  equations  of  this  report  can  be  used  to  study  the  complete  problem  for  both  symmetric 
and  asymmetric  masses.  Table  B-1  gives  the  necessary  parameters  for  the  projectile  studied  by 
Murphy.^  In  Figure  B-1,  we  plot  vs.  6)3/^,^  for  different  beam  damping,  =  0, 0.1,  0.2 . 

We  see  that  dynamic  instability  does  occur  near  when  ks  is  0.2.  Figure  B-2  shows  the 

damping  of  the  slow  frequency,  ^  ^  function  of  6)3 .  We  see  that  this  mode  is  imstable 
near  =i>2R  as  was  predicted  by  Reed.  In  Figure  B-3,  we  plot  and  vs.  . 
Resonance  occurs  for  p  =  \ .01  d)^,  and  this  large  beam  deflection  was  predicted  by  Reed.  This 

analysis  can  be  extended  to  other  beam  types,  pinned,  fixed-fixed,  etc.,  by  using  the  beam 
influence  coefficients  given  by  Chadwick  and  Murphy. 


*  Chadwick,  W.  R.  “Stability  of  Spinning  Shell.”  Unpublished  NSWC  Report,  Naval  Surface  Weapons  Center, 
Dahlgren,  VA,  September  1975. 

^  Murphy,  C.  H.  “Spinning  Projectile  Instability  Induced  by  an  Internal  Mass  Mounted  on  an  Elastic  Beam.”  ARL- 
MR-270,  U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  November  1995;  see  also  AIAA 
paper  92-4493CP,  August  1992. 

^  Reed,  H.  L.  “On  the  Motion  of  a  Shell  With  an  Internally  Mounted  Cantilever  Beam.”  BRL-TR-1039,  U.S.  Army 
Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  MD,  November  1957. 

^  Miuphy,  C.  H.  “Influence  of  Moving  Interior  Parts  on  the  Angular  Motion  of  Spinning  Projectiles.”  Journal  of 
Guidance  and  Control,  vol.  1,  no.  2,  pp.  117-122,  March-April  1978;  see  also  BRL-MR-2731,  U.S.  Army 
Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  MD,  February  1977. 
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Table  B-1.  Parameters  for  Spinning  Projectile  With  Interior  Beam-Mounted  Mass 


Intentionally  left  blank. 
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List  of  Symbols 


d 

Ej 


(EI)i 


gi 

Hj 

Ix 

Ixj 


It 


Iti 

m 


nij 


P 

Pj 

Q 

fm 

s 

V 


Xj 


projectile  diameter 

dimensionless  complex  lateral  location  of  cm  of  body  j 
with  respect  to  projectile  cm 

stiffness  of  beam  attached  to  body  j 

pV^S/2 

angular  momentum  vector  of  body  j 

axial  moment  of  inertia  of  projectile 
axial  moment  of  inertia  of  body  j 
transverse  moment  of  inertia  of  projectile 
transverse  moment  of  inertia  of  body  j 
dimensionless  length  of  body  j,  i.e.,  fineness  ratio 
projectile  mass 
mass  of  body  j 
Pi ,  projectile  spin 

axial  component  of  angular  velocity  of  body  j 
6  +  iif/ ,  complex  angular  velocity  of  body  1 
»22//m  ,  mass  ratio 
Tcd^jA 

magnitude  of  projectile  velocity 

dimensionless  axial  location  of  cm  of  body  j 
with  respect  to  the  cm  of  the  projectile 
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p 

k 

K 

e 

p 

<7 

G) 

&2 

(Of, 

(0„ 

¥ 

(4-4 ’4) 

(/V  /^  \ 


angle  of  attack  at  the  projectile  cm 
angle  of  sideslip  at  the  projectile  cm 
roll  angle  of  body  1 
frequency  of  k-th  mode 
damping  of  k-th  mode 
pitch  angle  of  body  1 
air  density 

common  value  of  beam  frequency  when  g)^  - 

^3{El\/[m,m2b^) 

K-th  frequency  for  free-free  uniform  elastic  beam 
,  rigid  projectile  zero-spin  frequency 
P  +  ia ,  complex  angle  of  attack  of  projectile 
K-th  eigenfunction  for  free-free  uniform  elastic  beam 
yaw  angle  of  body  1 
aerodynamic  force  exerted  on  body  j 

elastic  beam  force  exerted  on  body  j 
aerodynamic  moment  exerted  on  body  j 
elastic  beam  moment  exerted  on  body  j 
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(v^^. ,  Vyj ,  V^J  )  velocity  of  body  j 


{Xj ,  Yj ,  Zj )  location  of  cm  of  body  j  with  respect  to 
the  cm  of  body  1 

(rj,y ,  T^j )  pitch  and  yaw  angles  of  body  j  with  respect  to 
body  1 

r{z\  real  part  of  z 

l{z]  imaginary  part  of  z 

NOTES: 

-Tilde  superscript  denotes  quantity  related  to  elastic  beam. 
-Circumflex  superscript  denotes  aerodynamic  quantity. 

-e  subscript  denotes  earth-fixed  axes. 

-B  subscript  denotes  parameter  for  bent  projectile. 

-R  subscript  denotes  parameter  for  rigid  symmetric  projectile. 
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